Quasiparticles and the quantum fluctuations of local observables 
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The role that quasiparticles play in a strong interaction 
system with spontaneous symmetry breaking is examined. We 
find, using a non- perturbative cluster decomposition method, 
that the quasiparticles do not saturate the physical local ob- 
servables at small distances. The fermion number density 
serves as a clearcut example. A component due to localized 
random quantum fluctuations of the order parameter(s) in the 
vacuum state and the contributions of "quasiparticles" corre- 
sponding to other local minima of the effective potential is 
needed. At large distances, the ordinary quasiparticle picture 
emerges in the response of the system to classical background 
fields but the above mentioned component acts as a source 
for them. 
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In the self-consistent mean field and Hartree-Fock ap- 
proximations, the behavior of a quasiparticle is similar to 
an ordinary particle. A particle in a free or weakly inter- 
acting theory has definite mass; observables of the theory 
are saturated by the contributions of these particles. As 
a consequence, the fermion number density of a system 
with massive fermions is non-zero only after the chemical 
potential p exceeds the mass of the lightest fermions in 
the system. It is expected that the quasiparticles do the 
same at the mean field level. 

Quasiparticles are however not particles. The essen- 
tial difference between the mass of a particle and that 
of a quasiparticle lies in the fact that the mass of a par- 
ticle is a macroscopic parameter of the theory, with no 
quantum fluctuations, and the "mass" of the quasiparti- 
cle, being proportional to the order parameter, is a local 
variable. Such a local variable contains quantum fluc- 
tuations that do not decouple in the thermodynamical 
limit. As a result the fermion number density is not sat- 
urated solely by the contributions of the quasiparticles 
but contains a component due to local quantum fluctua- 
tions. This component, being non-propagating and dis- 
correlated with each other in different spacetime region, 
is called the dark component. 

Recent lattice studies of the spontaneous chiral sym- 
metry breaking in QCD at finite density found the emer- 
gence of baryon number at chemical potential below the 



expected threshold for particle production 
the non-quenched approach [Q. The origin of this effect 
is not yet well understood. It is found here that the dark 
component is responsible for such a behavior. 

For simplicity, the two flavor half bosonized NJL model 
is used for our investigation. The model Lagrangian den- 
sity with a chemical potential p is 



- p^/° — a — ij 5 f ■ ffj if) 



— (- 2 
2G 1 



(1) 



where Go is the coupling constant, a and ff are auxiliary 
fields and f is a set of three Pauli matrices. It has a 
SU(2) L x SU(2) R chiral symmetry. 

After performing the integration over the fermion fields 
if) and if), the generating functional W[J] can be written 
as 
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where "J" and "/" represent, collectively, the external 
fields and auxiliary fields respectively. The effective ac- 
tion S e ff is given by 
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where Sp denotes the functional trace. The operator 
iSp l [(j, ff, p] = ift + p/y° — a — vf'T ■ if is the inversed 
propagator of the fermions in the background auxiliary 
fields <j{x) and ff(a;). 

Let us define p[a, ir , p; x] = Tr7°(a; |Sir[<7, ff, x), 
where the trace "Tr" is over the internal degrees of free- 
dom of the fermions. The fermion number density of the 
model in a state | iff) determined by the initial and final 
configurations of the auxiliary fields is 
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where Z = J D [tj, ff] e 
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The Minkowski space- 
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time is not suitable to study the properties of the vacuum 
state using Eq. |] since the initial and final auxiliary field 
configurations are not specified. For a given set of con- 
figurations, the contributing intermediate states are not 
merely the vacuum states but are contaminated by the 
excited states. The usual procedure to project out the 
contributions of the vacuum state is go to the Euclidean 
spacetime in which the vacuum state has lowest energy. 
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In the mean field approximation, the vacuum phase 
is determined by minimizing the effective potential 
V e ff(cr,n) = — S e ff[cr, 0, /Lt]/f2 with fl — > oo the spacetime 
volume of the system, a spacetime independent and if as- 
sumed zero. Due to the chiral symmetry, the assumption 
if = does not loss generality. The phase of the sys- 
tem is determined by the condition 5V e ff(a,fi)/5a = 0. 
The solution for a is non-zero after the coupling constant 
Go is greater than a critical value Gq c . A non- vanishing 
a generates an effective mass for the fermions to form 
quasiparticles. 

The mean field fermion number density for the vac- 
uum is obtained from Eq. ^ by ignoring the functional 
integration over a and if and let a = a. The result is |J 
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which is non-zero only when (i > a just like the fermion 
number density of free massive particles with mass m = 
a. Therefore, the quasiparticle contributions saturate the 
fermion number density in the mean field approximation. 

The contributions of quantum fluctuations around the 
mean field tr are formally included in Eq. ||. The re- 
sults are commonly expressed as loop corrections to the 
fermion number density vertex. It is not attempted here. 

Instead of performing a loop expansion computa- 
tion of the fermion number density, we evaluate non- 
perturbatively the effects of the quantum fluctuations by 
"doing" the path integral. 

To proceed, the system under consideration is first 
putted in a Euclidean spacetime box of length L with 
periodic boundary conditions at its surfaces. The ther- 
modynamical limit is defined as the limit of L — > oo. In 
the thermodynamical limit, the extremal configuration 
dominates the path integral among those configurations 
of a and if that give divergent action in the thermody- 
namical limit. The contributing finite action quantum 
fluctuation configurations, the number of which is pro- 
portional to the spacetime volume Q — L 4 , are further 
classified into two categories: 1) correlated localized con- 
figurations, which are defined as the ones that approaches 
to the mean field configurations in the spacetime infinity 
and 2) correlated extended configurations, which are the 
ones that remain different from the mean field configura- 
tions at the spacetime infinity. 

For a given system, whether or not a configuration is 
an correlated extended configurations or is of localized 
onesQ is determined by dynamics. The on-shell ampli- 
tudes, are solutions of the "classical equation of motion" 
in the Euclidean spacetime 
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with S^ff the Euclidean effective action, / representing a 
or if fields. The superscript E shall be suppressed in the 
following. The set of the extended solutions to Eq. ^| are 
correlated ones. Albeit there are plenty of extended on- 
shell amplitudes in the Minkowski spacetime, there is no 
known one in the Euclidean one. We shall assume the ab- 
sence of them. The degree of correlation of an arbitrary 
extended configuration at different spacetime points is 
determined by the degree of their deviation from the ex- 
tended on-shell amplitudes for propagating excitations of 
the system. 

The correlation between the off-shell configurations at 
two different space time points decreases exponentially. 
For example, consider the field-field correlations or prop- 
agators 
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with T denoting time ordering. The off-shellness of these 
configurations is measured by their mass m. The cor- 
relation of these configurations at two different location 
X\ and X2 (in the Euclidean spacetime) decreases as fast 
as exp(— \x± — X2\m). So they can be decomposed into 
superposition of localized ones with a size of order 1/m. 
For these set of configurations, one can divide the space- 
time into cells with a dimension sufficiently larger than 
their correlation length. Then the contribution of this 
set of field configurations to the partition functional Eq. 
can be cluster decomposed to 
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with w k [J] = In z k [J] , z k [J] and w k [J] the corresponding 
partition functional of the kth cell. The full partition of 
the kth cell is 
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where S^f* [/] is the effective action of the kth cell and 
the path integration of / is over those ones that equals 
to the mean field value / outside of the kth cell but with 
arbitrary amplitudes inside the finite volume. For a given 
theory, instead of arbitrary division of the spacetime, it 
is expected that there is an optimal one with minimum 
volume n7 for each cell and yet has an error below a pre- 
determined one. We shall assume that such an optimal 
division of the spacetime into cells has already been found 
in the following discussion. 

Let us evaluate the contributions of the uncorrelated 
localized quantum fluctuations of the a and if fields to 
the vacuum fermion number density using Eqs. f§ | and 
|9| in the Euclidean spacetime. 

In the phase where the chiral SU(2)l X SU(2)r sym- 
metry is spontaneously broken down, the "chiral angle" 
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variable represented by 7? in the phase where (0 |tt| 0) = 
becomes massless following the Goldstone theorem. The 
correlation length of the 7? field becomes divergent in the 
chiral symmetrical limit. Therefore the field configura- 
tions of the Goldstone boson degrees of freedom contain 
the dominating on-shcll components that can be included 
by doing a loop expansion as usual. Such a loop ex- 
pansion contains no infrared divergences. Because the 
fermion number density p(a, ff, p; x) under study is chi- 
ral symmetric, which means that it does not depend on a 
spacetime independent global "chiral angle" ; it depends 
only on the derivatives of the "chiral angle" variables. 
The absence of a dependence of the fermion number den- 
sity on a global "chiral angle" guarantees the absence of 
the infrared divergences in the quantum corrections due 
to the Goldstone bosons. The configurations of the "chi- 
ral angle" , being on the edge of their shell and extended 
configurations in nature, are uniform and infinitesimal 
in amplitudes since it contains no infrared divergences 
and has a number of distinct modes proportional to the 
volume f2 of the system. They can not modify the quali- 
tative features of the quasiparticles. So, the tt variable in 
the vacuum fermion number density can be eliminated. 
It is treated as zero in the following discussions. 

The quantum fluctuation in the "mass" term, namely 
the chiral radius or order parameter represented by a 
(when (0 |7r| 0} = 0) has different characteristics due to 
the fact that it contains no on-shcll Euclidean configu- 
rations. These off-shell configurations have only short 
range correlations in spacetime. 

If the spacetime is divided into cells with their di- 
mension optimally determined, then the path integration 
within each cell can be done independently. This gives 
us a cluster decomposed partition functional of the form 
given by Eq. |^ with the partition functional for each cell 
computed independently of each other. 

The cluster decomposition property of the partition 
functional of the system reduces the full fermion number 
density of the vacuum given by Euclidean form of Eq. || 
to 
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where the cell labeled by k is the one that contains the 
spacetime point x, J Dk[a'] denotes integration over field 
configurations that approach the mean field value outside 
the cell and is the Euclidean effective action of the 
cell. 

Eq. [l(] is still too complicated to evaluate analytically. 
We make a further simplification by assuming that the 
functional integration of a' within a spacetime cell can 
be replaced by an ordinary integration over the spacetime 
averaged value of a' within that cell. It can be achieved 
by keeping the average of a' fixed while integrate over 
the rest degrees of freedom. The non-trivial part of it 
is in the assumption that after eliminating the rest of 



the degrees of freedom, the resulting effective potential 
V e ff remains, at least in form, the same as the original 
one. This procedure is in the same spirit as the renor- 
malization group analysis. In this way, Eq. [l^ reduces 
to 

1 f°° 
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where 5a =< a' > is the spacetime average of a' within 
the cell and the same reduction is also made to z[0]. The 
finiteness of w result in different qualitative behavior for 
p vac as a function of p. To explicitly see the difference, 
let us expand V e ff around ct, keeping only the leading 
quadratic term 

V eff (a + 5a) = V e ff (a) + A(a)Sa 2 + . . . (12) 

and using Eq. |s| for the trace term. If voA is sufficiently 
large, the result can be written as 
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It is clearly non-zero for any finite p below a since A is 

finite except in the limit of Go — > 0, which represents the 

free field case. The dependence of p\j a \ on p for a set of 

different values of a = wA is plotted in Fig. 1. Instead 
1/3 — 

of a sharp rise in p vac at p — a, p is non- vanishing all the 
way to p, = 0. This component of the fermion number 
density can certainly not be attributed to the contribu- 
tions of the quasiparticles. We find therefore that there 
is a dark component for the fermion number density that 
can not be accounted for by the quasiparticle contribu- 
tions. 

The NJL model has only one non-trivial vacuum phase 
since there is one independent absolute minimum in its 
effective potential. If the effective potential of the system 
contains a second local minimum with a higher energy 
density than the absolute minimum, then there are con- 
tributions from the local minimum to the fermion number 
density. The existence of such a contribution is another 
direct consequences of the cluster decomposability of the 
partition functional of the system. Models with a sec- 
ond minimum are studied in the literature. Those that 
have only one order parameter are represented by the 
Friedberg-Lee model B in which the two minima of the 
effective potential correspond to confinement and decon- 
finement phase of the model. Those that has two order 
parameters are introduced in Refs. |||| in which the (3 
and u> phases of the massless fermion system are studied. 

In the presence of a higher virtual phase that is sepa- 
rated from the actual phase of the system by a potential 
barrier, the fermion number density of the system is sat- 
urated by both the quasiparticles of the first phase and 
those ones of the second virtual phase together with their 
corresponding dark component discussed above. It has a 
form 
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where A = V}. 



V} f f > is the difference in energy 



'eff ~ V eff 

density between the virtual phase and the actual phase, 
p vaci and p va c 2 are the vacuum fermion number density 
of the actual vacuum phase and that of the virtual vac- 
uum phase respectively and w is the optimal volume of 
the spacetime cell between which the order parameters 
of the system are discorrelated. The contributions of the 
virtual phase to the fermion number density or to any 
local physical observables are non-perturbative effects. 
Attempts had been made in Refs. |||7]] to search for pos- 
sible other virtual phases J5|,|| of the strong interaction 
vacuum. 

The fermion number density discussed above are de- 
fined on a spacetime point. Such a precision is non- 
achievable in realistic observations. The physical observ- 
ables can be represented by a coarse-grained averaging of 
the form p vac = AN vac / An with An the smallest volume 
in spacetime that the observation apparatus can resolve 
and A7V„ ac the average number of fcrmions coherently 
produced by an external classical field K within that vol- 
ume. AN vac is not in general identical to d 4 xp vac 
when An >> vo. It is obtained from the partition func- 
tional W[J,p] by adding the external field K to p. K 
has a constant non-zero value only within spacetime vol- 
ume An, then AN vac = dW[0, p + K]/diK\ K=0 . In 
case when A57 << w, namely the precision of the ob- 
servation is much higher than the correlation length of 
the order parameter, the observed fermion number den- 
sity p V ac behaves in the same way as p vac - On the other 
hand, if An » zu, then the smallest spacetime cell that 
contributes to AN vac is a region of volume An rather 
than w. In this case, p vac is obtained from p vac by sub- 
stituting An for w. The effects of the dark component 
in the observed fermion number density are reduced as 
a result. When the resolution is sufficiently low, which 
means ASIA >> 1, the quasiparticle picture with only 
one vacuum , namely the actual one, reemerges in the 
response of the system to K . Therefore the dark compo- 
nent of the fermion number density is of transient nature. 

On the other hand, the K field radiated by the vac- 
uum fermion number density is of the form K(x) — 
J d 4 x'G(x,x')p va c(x') with G(x,x') the Green function 
for K(x). It is p vac that is the source for K(x) rather 
than p vac no matter how slow the resulting K(x) varies 
in spacetime. Thus the effects of the dark component are 
indirectly observable even in low energy processes. 

Albeit the fermion number density is discussed, the 
results of the letter are also applicable to other local ob- 
servables like the energy density. The finding may have 
implications on the dark matter problem in cosmology 
since it implies that the apparent matterless space at the 
macroscopic level may reveal itself of matter effects in 
low energy gravitational processes even when p is below 



the baryonic particle production threshold and the p = 
energy density subtracted. This is because gravitational 
fields couple locally to the source matter fields that con- 
tain the localized dark component. 

In summary, it is found that the quasiparticle picture 
is insufficient in describing quantum processes of strong 
interaction at both small and large distance scales. The 
random quantum fluctuations of order parameters of the 
system can change the behavior of the system qualita- 
tively. The finding may be relevant to our understanding 
of the dark matter problem in astronomy. Other impli- 
cations in hadron physics are also worth exploring in the 
future. 
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FIGURES 



Fig. 1 The dependence of p vac on the chemical potential 
p. Here a = wA. The unit for the dimensional quantities 
are GeV. Solid lines represent the case of free theory with 
mass and 0.5 respectively. Other lines represent the 
results for the massless NJL model with different strength 
of local fluctuations characterized by the a of the order 
parameter. 
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